ABSTRACT 

New  human  estimates  for  chlorine 
inhalation  lethality  as  a  function  of 
exposure  duration  were  derived 
via  a  review  and  statistical  analysis  of 
exisfing  mammalian  lefhalify  data.  Such 
estimates  are  needed  to  support  risk  assess¬ 
ments  and  casualty  predictions  involving 
airborne  releases  of  chlorine.  At  present, 
casualty  predictions  for  such  releases  are  af 
odds  with  what  has  been  observed  histor¬ 
ically;  the  predicted  downwind  hazard  area 
has  often  been  much  larger  than  what  was 
actually  observed.  Either  the  present  esti¬ 
mates  for  median  lefhal  dosages  (LCT50)  are 
foo  low,  fhe  currenfly  popular  atmospheric 
transport  and  dispersion  (ATD)  models 
cannot  adequately  model  chlorine  releases, 
or  both. 

LCT50  and  quantal  response  data  were 
analyzed  for  eight  species  (mouse,  rat, 
guinea  pig,  rabbit,  cat,  dog,  goat  and  sheep), 
for  exposure  durations  from  8  fo  235  min- 
ufes.  The  base  10  probit  slope  (concentra¬ 
tion)  was  estimated  via  the  weighted 
average  of  experimenfally  measured  slopes 
in  mammalian  lethality  studies.  Resulting 
human  lethality  (military)  estimates  as 
a  function  of  exposure  duration  were 
expressed  via  the  toxic  load  model.  General 
population  estimates  were  derived  from  the 
military  estimates  using  the  mathematical 
method  of  Crosier  (2007). 

Previous  human  esfimafes  were  re¬ 
viewed  and  compared  to  the  new  general 
population  estimate.  The  impact  of  fhe 
new  estimate  was  evaluated  through  a  se¬ 
ries  of  fransport  and  dispersion  modeling 
runs  for  fhe  catasfrophic  accidental  release 
of  50  tons  of  chlorine  from  a  tanker  car. 
The  sensitivity  of  downwind  hazard  dis- 
fances  was  also  invesfigafed  as  a  function 
of  median  lethal  toxic  load  (TL),  toxic  load 
exponent  and  probit  slope  values. 


INTRODUCTION 

Chlorine  is  a  greenish-yellow,  highly  re¬ 
active  halogen  gas  that  has  a  pungent,  suf¬ 
focating  odor.  Because  of  its  extensive 
industrial  use  and  toxicity,  it  has  been  the 
subject  of  many  toxicological  studies,  and 
a  number  of  human  toxicity  estimates  have 
been  derived.  Several  previous  efforts  cen¬ 
tered  on  toxicity  estimates  for  protecting  a 


population  from  accidental  releases.  How¬ 
ever,  many  risk  assessment  applications  re¬ 
quire  an  objective  estimate  of  the  actual 
toxicity  and  not  a  conservative,  protective 
estimate.  At  present,  casualty  predictions 
for  accidenfal  chlorine  releases  are  af  odds 
wifh  what  has  been  observed  historically. 
Either  the  present  toxicity  estimates  are 
too  low  (i.e.  potency  is  over  estimated),  the 
currently  popular  atmospheric  transport 
and  dispersion  (ATD)  models  cannot  ade¬ 
quately  model  chlorine  releases,  or  both. 

The  purpose  of  this  study  is  to  reevalu¬ 
ate  the  human  lethality  estimate  as  a  fimction 
of  duration  for  chlorine  inhalation  (IH)  ex¬ 
posures  involving  a  healfhy  (military-fype) 
human  populafion.  An  allomefric  modeling 
approach  was  used  to  extrapolate  a  human 
estimate  from  applicable  mammalian  foxic- 
ity  data.  The  estimate  also  provides  a  probit 
slope  for  a  healthy  adult  population.  Gen¬ 
eral  population  estimates  were  developed 
from  the  military  estimates  using  Crosier 
(2007). 

STATISTICAL  BACKGROUND 
Allometric  Modeling 

The  dependence  of  a  biological  variable 
X  on  body  mass  M  is  typically  characterized 
by  an  allometric  scaling  law  of  fhe  form: 

X  =  XoM^  [1] 

where  b  is  fhe  scaling  exponent  and  Xg  a  fit¬ 
ted  constant  (West  et  al.,  1997).  These  two 
coefficients  can  be  estimated  via  the  linear 
regression  fit  of  X  versus  M  on  a  log-log 
scale,  where  fhe  slope  and  y-intercept  of 
the  resulting  line  equal  b  and  log(Xo),  re¬ 
spectively. 

In  inhalation  toxicology,  the  dose  re¬ 
ceived  by  an  animal  is  directly  related  to 
the  amount  of  air  inhaled  (Phalen,  1984, 
Salem,  1987,  Snipes,  1988).  Minute  volume 
(Vm)  is  a  commonly  used  respiratory  mea¬ 
surement,  which  equals  the  volume  of  air 
inspired  per  minute.  A  relatively  recent  re¬ 
view  (and  analysis)  of  mammalian  Vm  dafa 
(and  species  body  mass,  M)  is  available  from 
Bide  et  al.  (2000).  [Table  1  lists  some  of  fheir 
average  Vm  and  M  values.]  Based  upon  allo¬ 
metric  scaling  law,  Vm  should  scale  to  the 
0.75  power,  but  they  found  thaf  Vm  scales 
to  the  0.809  power  for  non-anesfhetized, 
adult  mammals: 
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Table  1.  Summary  of  Minute  Volumes  (Non-Anaesthetized  Animals)’. 

Species 

Mass 

(kilograms) 

Minute  Volume 
Expt.  Average 
(liters) 

Minute  Volume 
Allometric  Model  Fit 
(liters) 

Vm  Ratio 
(Expt  /  Fit) 

Mouse 

0.0207 

0.0310 

0.0217 

1.431 

Rat 

0.261 

0.203 

0.168 

1.206 

Guinea  Pig 

0.323 

0.182 

0.200 

0.910 

Rabbit 

2.77 

1.453 

1.138 

1.277 

Cat 

3.36 

0.845 

1.330 

0.635 

Dog 

13.0 

4.230 

3.975 

1.064 

Goat 

36.9 

9.600 

9.243 

1.039 

Sheep 

52.6 

13.90 

12.31 

1.129 

*  adapted  from  Bide  et  ah,  2000 


Vm  =  (0.499)  [2] 

where  is  in  liters/ min  and  M  is  in  kg.  Using 
Equation  [2],  they  calculated  a  of  15.5  liters/ 
min  for  70  kg  man  engaged  in  mild  activity — 
not  a  resting  volume.  This  Vm  is  very  close  to 
the  standard  Vj^  of  15  liters/min  used  by  DoD 
as  a  basis  for  their  toxicity  estimates  for  a  70 
kg  soldier  (US  Army  Chemical  School,  2005). 

For  typical  whole  body  exposures,  the  exact 
inspired  amount  is  usually  not  measured.  In  its 
place,  a  nominal  dose  (mg)  can  be  calculated  by 
multiplying  the  lethal  inhalation  dosage,  (mg- 
min/m^)  (estimated  from  experimental  toxicol¬ 
ogy  data)  by  the  minute  volume  (m^/min): 

LD50  =  Vm  X  LCT50  [3] 

where  LD50  is  the  nominal  dose  lethal  to  50% 
of  exposed  individuals  and  LCT50  is  the  lethal 
concentration-time  (or  inhalation  dosage)  to  50% 
of  exposed  individuals. 

Toxicity  is  a  function  of  many  factors,  and 
some  of  these  {i.e.  metabolic,  circulatory,  etc.) 
may  also  scale  allometrically  in  their  effect.  Al- 
lometric  modeling  can  be  used  to  empirically 
account  for  the  cumulative  effect  of  these  factors 
and  scale  the  nominal  LD50  (mg)  as  a  function  of 
species  body  mass,  as  shown  in  Equation  [4].  A 
human  estimate  can  be  obtained  from  mamma¬ 
lian  data  by  extrapolating  to  a  70  kg  human. 

LD5o  =  XoM'’ 

or  [4] 

log(LD5o)  =  log(Xo)  -t  fclog(M) 


In  Equation  [4],  both  the  regressor  term  (LD50) 
and  the  predictor  term  (M)  are  highly  correlated 
with  the  parameter,  Vm-  Controversy,  dating 
back  to  1897  (Pearson,  1897),  exists  over  regres¬ 
sions  involving  predictor  and  regressor  vari¬ 
ables  sharing  a  common  term,  and  continues 
to  this  day  (Albrecht  et  al.,  1993;  Bradshaw  & 
Radbill,  1987;  Dunlap  et  al.,  1997;  Firebaugh, 
1986,  1988;  Firebaugh  &  Gibbs,  1985,  1987; 
Kraft,  1987;  Kuh  &  Meyer,  1955;  Macmillan  & 
Daft,  1980;  Prairie  &  Bird,  1989;  Smith,  2005; 
Vanderburgh,  1998;  Yule,  1910) — it  is  thought 
that  such  regressions  return  spurious  correla¬ 
tions.  The  statistical  community  has  rejected 
the  notion  that  these  t5rpes  of  regression  are  sta¬ 
tistically  invalid,  but  many  individuals  in  the 
greater  scientific  community  have  not  kept 
abreast  of  these  statistical  advances  (Firebaugh 
&  Gibbs,  1985, 1987;  Prairie  &  Bird,  1989). 

Inhalation  Dose-Response  Statistics 
and  the  Toxic  Load 

The  statistical  properties  of  dose-response 
distribution  have  been  extensively  reviewed 
(Finney,  1971;  Neumann  &  Kimmel,  1998) 
and  summarized  (Salem,  1987;  Crosier  & 
Sommerville,  2002).  Inhalation  dosage  tradi¬ 
tionally  has  been  expressed  as  the  product  of 
vapor  concentration  (C)  and  duration  (T),  or 
CT.  The  distribution  of  effective  CTs  for  a  ho¬ 
mogeneous  population  is  usually  lognormal 
and  described  by: 

Z  =  m  [log(LCTxx)  -  log(LCT5o)]  [5] 
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where  Z  is  the  standard  normal  random  vari¬ 
able,  m  is  the  probit  (or  Bliss)  slope  (and  equal 
to  inverse  of  the  standard  deviation),  LCTxx 
and  LCT50  are  the  CTs  that  will  produce  lethality 
in  XX%  and  50%,  respectively,  of  the  exposed  in¬ 
dividuals.  Z  is  a  function  of  the  cumulative 
probability,  with  the  16,  50  and  84%  response 
level  corresponding  to  Z  equal  to  -1,  0  and  1,  re- 
spechvely  (Box  et  ah,  1978;  Hine  &  Montgomery, 
1990;  MINITAB®  Statistical  Software,  2006). 
Toxicologists  traditionally  use  base  10  logarithms 
to  calculate  the  probit  (Bliss)  slope  (Bliss,  1934; 
Finney,  1971),  while  engineers  often  use  natural 
logarithms  (Fairhurst  &  Turner,  1993;  Ferguson 
&  Hendershot,  2000;  The  Major  Hazards  Assess¬ 
ment  Panel,  1988;  Yee,  1996a). 

When  the  LCT50  is  constant  with  respect  to 
T,  it  is  said  to  obey  Haber's  Rule  (Haber,  1924). 
However,  Haber's  Rule  is  more  of  an  exception 
than  a  true  rule.  Thus,  a  new  term.  Toxic  Load 
(TL),  has  been  developed  and  extensively  used 
in  inhalation  toxicology  to  account  for  time- 
dependent  toxicity  (Mannan,  2005;  Ride,  1995; 
Sommerville  et  ah,  2005;  Yee,  1996b;  Yee  &  Ye, 
1996).'  TL  is  normally  expressed  as  C”T,  with  n 
being  the  toxic  load  exponent  (TLE).  LTL50  re¬ 
fers  to  the  TL  that  will  produce  lethality  in 
50%  of  exposed  individuals.  When  n  =  1,  the 
chemical  follows  Haber's  Rule  (Sommerville 
et  al.,  2005).  For  acute  inhalation  of  irritant  gases, 
n  tends  to  be  greater  than  unity  and  is  often  of 
the  order  of  two  in  value  (Mannan,  2005).  The 
same  statistical  theory  used  for  dosages  (Finney, 
1971)  has  also  been  adopted  for  use  with  toxic 
load  (Fairhurst  &  Turner,  1993;  Griffiths,  1991; 
Schubach,  1995,  1997;  Withers  &  Lees,  1985a,b; 
Yee,  1996a,b). 

Although  the  normal  distribution  is  contin¬ 
uous,  quantal  data  (response  versus  no  response) 
are  used  to  fit  Equation  [5]  (Eirmey,  1971, 1978), 
through  the  use  of  probit  analysis  and  maxi¬ 
mum  likelihood  estimation  (MLE)  (Eairhurst  & 
Turner,  1993;  Einney,  1971;  Eox,  1997).  The  fol¬ 
lowing  equation  is  fitted  via  probit  analysis/ 
MLE  for  vapor  toxicity  studies  (Eairhurst  & 
Turner,  1993;  Eirmey,  1971;  Yee,  1996a,b): 

Yn  =  (Yp  -5)=ko+  fcclogC  +  krlogT  [6] 

where  is  a  normal  probability  unit  (or  nor- 
mit),  Yp  is  a  probability  unit  (or  probit),  the  k’s 
are  fitted  coefficients,  C  is  vapor  concentration. 


and  T  is  exposure  time.  The  constants  kc  and 
kp  are  the  probit  slopes  for  concentration  and 
time,  respectively.  The  TLE  (n)  equals  the  ratio 
of  these  two  slopes,  or  {kc  /  kp).  Some  studies 
report  a  probit  slope  for  the  TL  (Eairhurst  & 
Turner,  1993;  Eranks  et  al.,  1996;  Marshall, 
1989;  Poblete  &  Lees,  1984;  ten  Berge  &  van 
Heemst,  1983;  Withers  &  Lees,  1985a,b). 

When  fitting  Equation  [6],  all  variability  in 
the  data  will  contribute  to  the  estimate  for  m, 
be  it  from  variance  due  to  individual  suscepti¬ 
bilities  (the  true  goal),  batch  effects,  experimen¬ 
tal  error,  differences  in  durations,  compilation 
from  many  sources  (Eranks  et  al.,  1996),  etc. 
The  heterogeneity  introduced  by  outside  sour¬ 
ces  of  variance  will  artificially  lower  the  probit 
slope.  The  precision  of  slope  estimation  typi¬ 
cally  improves  as  more  experimental  subjects 
are  used.  Eor  studies  involving  more  than  one 
exposure  duration,  fitting  Equation  [6]  with  T 
treated  as  a  factor  (instead  of  as  a  covariate)  will 
produce  the  best  overall  probit  slope  estimate 
for  the  whole  study  (Benton  et  al.,  2007;  Hulet 
et  al.,  2006).  Equation  [6]  can  also  be  written  in 
the  following  form: 

YN  =  K+kTplog{TL)  [7] 

with  kg  and  kpp  being  fitted  coefficients,  and  kpp 
being  the  probit  slope  on  a  toxic  load  basis.  Also, 
kpi  equals  the  probit  slope  for  time,  kp. 

The  TL  relationship  is  empirical,  rather  than 
theoretical  (Eairhurst  &  Turner,  1993;  Griffiths, 
1991;  Yee,  1996a,b).  Thus,  Equations  [6]  and  [7] 
needs  to  be  empirically  derived  on  an  individ¬ 
ual  toxicant  basis  from  experiments  where  both 
C  and  T  are  varied  (ten  Berge  et  al.,  1986).  The 
value  for  n  may  offer  some  insight  into  the  influ¬ 
ence  of  underlying  toxicological  mechanisms 
(Griffiths,  1991).  Eor  instance,  it  can  be  argued 
that  if  n  >  1,  then  significant  detoxification 
may  be  occurring  during  the  exposure,  since 
the  required  effective  dosage  is  increasing  as 
the  dosage  is  incurred  over  a  longer  duration. 

When  performing  an  allometric  regression 
fit  of  inhalation  dose  data.  Equation  [4]  is  ex¬ 
panded  to  the  following  form  to  account  for 
the  time  dependency  of  the  toxicity: 

log(LD5o)  =  log(Yo)  +  fcMlog(M)  +  kplogiT) 

[8] 
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Using  Equations  [3]  and  [8],  it  can  be  shown  al¬ 
gebraically  that  the  TLE  (n)  equals  [1  /  (1  -  kj-)]. 
After  the  appropriate  substitutions.  Equation 
[8]  can  be  written  in  a  TL  form  (Equation  [7]). 

Statistical  Treatment  of 
Historical  Data 

Proper  evaluation  of  historical  data  is  criti¬ 
cal.  The  studies  were  often  done  in  different 
decades,  by  different  researchers  in  different 
laboratories.  With  linear  regression  analysis, 
the  use  of  weights  technically  permits  assigning 
a  numerical  value  to  the  quality  of  a  datum,  but 
how  does  the  analyst  properly  estimate  the 
values  for  weighting?  A  solution  to  this  problem 
is  the  use  of  one  of  many  robust  regression 
methods  that  have  been  developed  (Box  & 
Draper,  1987;  Eox,  1997;  Myers  &  Montgomery, 
1995). 

Eor  the  present  study,  the  iteratively  re¬ 
weighted  least  squares  (IRES)  method  (Box  & 
Draper,  1987;  Eox,  1997;  Myers  &  Montgomery, 
1995)  was  used.  The  IRES  method  does  not 
involve  making  an  a  priori  decision  on  what 
weight  values  to  use.  Instead,  through  an  itera¬ 
tive  process,  the  weight  of  a  datum  is  calculated 
as  some  function  of  its  residual  from  the  previ¬ 
ous  iterative  step.  Often,  the  standardized  resid¬ 
ual  is  used  (Box  &  Draper,  1987),  but  for  this 
study  the  studentized  deleted  residual  was  used 
instead.  The  studentized  deleted  residual  of  an 
observation  uses  a  predicted  value  from  a  re¬ 
gression  fit  omitting  the  observation,  and  it  is 
more  sensitive  to  outliers  than  the  standardized 
residual.  Weights  based  on  deleted  residuals 
will  be  less  than  those  based  on  standardized 
residuals  (all  other  factors  being  equal). 

Military  Subpopulation  versus 
General  Population  in 
Toxicity  Sensitivity 

It  is  assumed  that  only  healthy  animals  were 
used  in  the  reviewed  studies.  Consequently, 
when  scaling  from  smaller  mammals  up  to 
humans,  the  derived  estimate  should  correspond 
to  a  healthy  human  subpopulation,  defined  in 
this  case  as  being  healthy  and  fit  enough  for  mil¬ 


itary  service.  However,  this  does  not  address  the 
issue  of  the  more  vulnerable  portions  of  the  gen¬ 
eral  population.  Several  studies  have  attempted 
to  account  for  such  sensitivity  differences  in  dif¬ 
ferent  populations  (Eisenberg  et  al,  1975;  Harris 
&  Moses,  1983;  MHAP,  1988;  NRC,  2004;  ten 
Berge  &  van  Heemst,  1983;  Withers  &  Lees, 
1985a,b).  Eor  this  study,  it  is  assumed  that  the 
sensitive  subpopulation  can  be  adequately  mod¬ 
eled  using  Equation  [5]  if  a  probit  slope  and  me¬ 
dian  effechve  dose  for  the  general  (or  whole) 
population  were  used.  The  vulnerable  subpopu¬ 
lation  would  be  accounted  for  by  using  the  lower 
tail  of  the  beU  distribution.  Thus,  the  key  is  to 
change  the  basis  (from  healthy  subpopulahon 
to  general  population)  of  the  probit  slope  and 
median  effechve  dose  estimates  that  were  de¬ 
rived  using  the  mammalian  data. 

In  a  series  of  reports  (Crosier  &  Sommerville, 
2002;  Crosier,  2003,  2007),  ECBC  has  developed 
a  mathematical  method  for  the  conversion  from 
one  population  basis  to  another,  with  Crosier 
(2007)  presenting  the  final  version  of  the 
method.  The  premise  of  the  approach  is  that 
mathematically  the  distribution  of  log(doses) 
for  a  healthy  subpopulation  is  located  com¬ 
pletely  within  the  distribution  formed  by  the 
general  populahon  {i.e.  the  smaller  bell  curve  is 
inside  the  larger  bell  curve).  A  major  advantage 
of  Crosier  (2007)  is  that  it  reduces  the  amount 
of  subjechvity  that  has  been  used  in  the  past  in 
accounting  for  populahon  differences. 

Military  Chemical  Warfare  Defense 
and  the  Toxic  Load  Model 

Use  of  the  Toxic  Load  Model  is  a  recent  de¬ 
velopment  in  military  toxicology,  being  offi¬ 
cially  adopted  by  the  Department  of  Defense 
(DoD)  in  2005  (US  Army  Chemical  School, 
2005).  The  demise  of  the  offensive  chemical  war¬ 
fare  (CW)  program  (and  the  resulting  change 
in  mission)  necessitated  such  a  course.  Using 
Haber's  Rule  was  sufficient  for  offensive  CW 
operations,  where  the  exposure  durations  of 
interest  were  from  2  to  10  minutes.  However, 
CW  defense  is  now  the  primary  focus.  Both 
short  and  long  exposure  durations  (from  a  few 
minutes  to  hours)  need  to  be  considered-for 
example: 
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•  Protective  equipment  tested  against  the 
higher  effective  concentrations  associated 
with  short  exposure  durations 

•  Decontamination  standards  (among  other 
items)  based  on  the  lower  effective  con¬ 
centrations  associated  with  longer  exposure 
durations 

Proper  risk  optimization  involving  levels  of 
CW  profection  requires  an  accurafe  knowledge 
of  the  TLE.  The  issue  is  illustrated  in  Figure  1, 
with  a  comparison  of  two  different  LCT50  time 
relationships:  Haber's  Rule  agent  {n  =  1)  versus 
a  TL  agent  having  an  n  >  1.  Decision  making 
using  Haber's  rule  for  a  TL  agent  (or  vice  versa) 
may  produce  an  inaccurate  assessment. 

HISTORICAL  CHLORINE  DATA 
AND  HUMAN  TOXICITY  ESTIMATES 

Toxic  Action  of  Chlorine  Inhalation 
on  Animals  and  Humans 

Chlorine  is  an  acute  respiratory  irritant  in 
both  animals  and  humans  and  produces  a  simi¬ 


lar  spectrum  of  foxic  responses  (from  minor  irri¬ 
tation  to  death),  in  all  species,  as  the  exposure 
concentration  increases  (MHAP,  1988;  Mannan, 
2005;  NRC,  2004;  Turner  &  Fairhurst,  1990). 
Chlorine  reaches  the  lungs  because  it  is  only 
moderately  soluble  in  water  and  it  is  not  totally 
absorbed  in  the  upper  respiratory  tract  at  high 
concentrations  (NRC,  2004).  It  has  been  con¬ 
cluded  by  the  Major  Hazards  Assessment  Panel 
(MHAP)  that  there  is  little  basis  for  supposing 
thaf  fhe  sensifivity  of  humans  fo  chlorine  is  sig¬ 
nificantly  different  from  that  of  small  mammals 
(MHAP,  1988;  Mannan,  2005). 

Previous  Human  Lethality  Estimates 
for  Chlorine  Inhalation  Toxicity 

There  are  several  excellenf  reviews  of  the 
body  of  chlorine  inhalation  toxicity  information 
available.  One  of  the  more  recent  is  Mannan 
(2005),  as  well  as  the  1987  monograph  by  the 
MHAP  sponsored  by  the  Institution  of  Chemi¬ 
cal  Engineers.  Human  eshmafes  from  five  previ¬ 
ous  studies  were  investigated  for  comparison 
with  the  present  work: 


U) 


Haber's  Rule 

—  O  Toxic  Load  Model  (with  n  >1) 


Value  for  toxic  load  exponent  (n) 
and  constant  (k„)  is  agent  and 
scenario  specific 


Haber's  Law 
Overestimates  Toxicity 


Exposure  Duration  (minutes) 


Figure  1.  Importance  of  Establishing  the  Time  Dependence  of  IH  Toxicity. 
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(1)  Eisenberg  et  al.  (1975) 

(2)  Perry  &  Articola  (1980) 

(3)  ten  Berge  &  van  Heemst  (1983) 

(4)  Rijnmond  Report  (1982)/Harris  &  Moses 

(1983) 

(5)  Withers  &  Lees  (1985a,b) 

All  five  studies  recognized  that  chlorine  inhala¬ 
tion  toxicity  varies  with  duration  and  used  the 
TL  model  to  describe  the  relationship  (see  "In¬ 
halation  Dose-Response  Statistics  and  the  Toxic 
Load").  Nonetheless,  these  studies  do  differ  in 
their  philosophies  on  how  to  best  extrapolate 
from  animal  data  to  a  human  toxicity  estimate. 
Although  some  (MHAP,  1987;  Withers  &  Lees, 
1985a)  have  observed  that  the  LC50  varied  as 
a  function  of  species  size,  there  are  no  known 
studies  that  have  actually  attempted  allometric 
modeling  of  the  existing  database. 

Eisenberg  et  al.  (1975)  and  Perry  &  Articola 
(1980)  are  related  works  that  were  written  in 
support  of  the  US  Coast  Guard's  Vulnerability 
Model-chlorine  was  one  of  many  chemicals 
for  which  human  toxicity  estimates  were  col¬ 
lected  or  derived.  These  human  estimates  lie 
well  below  the  LCT50S  of  the  underlying  exper¬ 
imental  animal  data,  and  several  authors  do  not 
recommend  their  estimates  for  risk  assessment 
(Harris  &  Moses,  1983;  MHAP,  1987;  ten  Berge 
&  van  Heemst,  1983;  Withers  &  Lees,  1985a,b). 

ten  Berge  &  van  Heemst  (1983)  assumed 
that  humans  can  be  modeled  as  being  as  sensi¬ 
tive  as  the  average  mammal,  which  would  place 
the  human  between  the  mouse  (most  sensitive) 
and  the  dog  (least  sensitive).  Their  estimate  of 
the  human  response  variability  was  for  a  general 
population.  The  toxicity  time  dependence  was 
based  upon  mouse  data  of  Bitron  &  Aharonson 
(1978)  and  human  volunteer  irritation  dose- 
response  relationship  reported  by  Anglen  (1981). 

Harris  &  Moses  (1983)  did  not  make  any  ex¬ 
plicit  statement  on  what  animal  model  they 
used  to  base  their  estimates.  Since  their  human 
estimate  curve  is  very  similar  to  that  of  ten 
Berge  &  van  Heemst  (1983)  it  may  be  inferred 
that  Harris  &  Moses  assumed  that  human  toxic¬ 
ity  can  be  equated  to  that  of  an  average  experi¬ 
mental  animal.  Harris  &  Moses  believed  that 
more  variability  exists  in  the  human  general 
population  response  (based  upon  the  historical 
record)  than  was  originally  estimated  by  the 


US  Coast  Guard  studies.  However,  the  Harris 
&  Moses  estimate  for  response  variability  is  less 
than  that  of  ten  Berge  &  van  Heemst.  Thus,  al¬ 
though  the  predicted  LGT50S  from  the  two  esti¬ 
mates  are  practically  identical  (for  all  exposure 
durations),  predicted  dosages  for  percent  effect 
levels  other  than  50%  will  increasingly  differ  as 
the  two  extremes  (zero  and  100%  response)  are 
approached.  Harris  &  Moses  and  the  Rijnmond 
Report  (1982)  have  identical  human  estimates 
and  are  treated  as  synonymous  by  Mannan 
(2005). 

Withers  &  Lees  (1985a,b)  gave  greatest 
weight  to  the  dog  data  of  Underhill  (1920)  in 
the  development  of  their  human  estimate,  but 
did  not  neglect  to  interpret  the  dog  data  in  light 
of  the  other  animal  data.  They  also  developed 
human  estimates  as  a  function  of  the  type  of  hu¬ 
man  population  (regular,  average  and  vulnera¬ 
ble)  and  two  levels  of  physical  activity  (base 
and  standard,  corresponding  to  minute  vol¬ 
umes  of  6  and  12  liters/ min,  respectively).  All 
of  the  Withers  &  Lees  estimates  are  lower  than 
those  for  the  general  population  estimates  from 
ten  Berge  &  van  Heemst  (1983)  and  Harris  & 
Moses  (1983),  with  the  difference  increasing  as 
the  duration  increases. 

HISTORICAL  RECORD  OF  CHLORINE 
RELEASES 

The  historical  record  of  chlorine  industrial 
accidents  has  been  well  documented  (Harris  & 
Moses,  1983;  MHAP,  1987;  Mannan,  2005; 
Marshall,  1977).  The  most  notable  statistic  for 
the  reported  incidents  of  accidental  outdoor  re¬ 
leases  (with  release  amounts  up  to  50  metric 
tons)  is  that  fatalities  (with  one  exception)  occur 
within  about  400  meters  of  the  actual  point  of  re¬ 
lease,  and  more  often  within  250  meters  (MHAP, 
1987;  Mannan,  2005).  The  largest  number  of  fa¬ 
talities  for  an  outdoor  chlorine  release  is  around 
60  and  resulted  from  a  ruptured  storage  tank  in 
Romania  in  1939  (Marshall,  1977).  It  has  been  es¬ 
timated  that  there  are  (on  average)  0.3  to  0.8  fa¬ 
talities  per  metric  ton  of  an  accidental  industrial 
release  (Harris  &  Moses,  1983;  MHAP,  1987; 
Marshall,  1977). 

Marshall  (1977)  made  the  anecdotal  obser¬ 
vation  that  no  civilian  casualties  were  reported 


Page  64 


Military  Operations  Research,  V15  N3  2010 


REVIEW  AND  ASSESSMENT  OF  CHLORINE  MAMMALIAN  LETHALITY  DATA  AND  THE 

DEVELOPMENT  OF  A  HUMAN  ESTIMATE 


from  the  first  large  scale  chemical  warfare  attack 
at  the  Second  Battle  of  Ypres  on  22  April  1915  in 
World  War  I  (Heller,  1984).  The  Germans  re¬ 
leased  168  tons  of  chlorine  over  a  four  mile 
front,  producing  an  estimated  5000  military  fa¬ 
talities.  The  distance  between  the  German  and 
Allied  lines  were  roughly  200  to  600  meters  in 
the  vicinity  of  the  attack.  The  Allied  lines  were 
located  about  four  miles  (or  6  km)  upwind  from 
the  town  of  Ypres,  which  was  still  occupied  by 
many  civilians.  It  might  be  inferred  from  the 
lack  of  evidence  that  the  chlorine  cloud  never 
reached  the  town,  since  it  is  likely  wartime  pro¬ 
paganda  would  have  made  the  most  of  any  ci¬ 
vilian  deaths  had  they  occurred. 

It  is  well  known  that  ATD  models,  in  gen¬ 
eral,  often  greatly  overestimate  the  downwind 
hazard  resulting  from  chemical  releases,  and 
modeling  of  chlorine  is  no  exception  (Marman, 
2005).  There  are  many  possible  factors  for  this: 
inadequacy  of  ATD  model  algorithms  (for  dense 
gas  releases,  terrain  effects,  etc.);  incomplete 
knowledge  of  the  atmospheric  chemistry  of  chlo¬ 
rine;  the  source  term  for  the  chlorine  gas  cloud; 
and  the  toxicity  estimates. 

Based  upon  the  above  historical  record,  any 
chlorine  toxicity  estimate  that  predicts  a  signifi¬ 
cant  number  of  fatalities  farther  than  a  kilometer 
downwind  from  the  release  point  should  be 
considered  suspect.  A  significant  number  of 
predicted  fatalities  between  500  and  1000  meters 
should  be  cause  for  concern.  A  case  in  point  is 
that  of  Franks  et  al.  (1996).  Using  a  toxicity  esti¬ 
mate  from  the  Health  and  Safety  Executive 
(UK)  (Turner  &  Fairhurst,  1990),  Franks  et  al.  pre¬ 
dicted  via  ATD  modeling  that  the  downwind 
range  of  a  LC95  would  be  from  400  to  800  meters 
for  a  chlorine  puff  release  of  1000  kg  (a  small 
amoimt  for  an  industrial-t5q)e  accident),  and  for 
the  LC50,  from  650  to  1100  meters. 


DATA  ACQUISITION 

An  extensive  literature  review  was  con¬ 
ducted  to  find  existing  LCT50  estimates  for  acute 
inhalation  exposures  to  chlorine.  Whenever  pos¬ 
sible,  the  underlying  raw  quantal  response 
data  were  also  obtained.  For  the  initial  identifi¬ 
cation  of  relevant  studies,  three  reference  works 
were  consulted:  the  Acute  Exposure  Guideline 


Levels  (AEGL)  for  chlorine  (NRG,  2004);  the 
Third  Edition  of  Lee's  Loss  Prevention  in  the 
Process  Industries  ("Toxic  Releases",  Chapter 
18)(Marman,  2005);  and  the  Chlorine  Toxicity 
Monograph  published  by  the  Major  Hazards 
Assessment  Panel  (MHAP,  1988).  The  search 
was  subsequently  expanded  beyond  this  point. 

A  total  of  29  individual  LCT50S  were  col¬ 
lected  (or  derived  from  raw  data)  from  18  stud¬ 
ies  and  sources  for  eight  species  [see  Table  2]. 
Whenever  possible,  original  sources  were  found 
for  each  LCT50.  In  some  cases,  data  were  discov¬ 
ered  in  previously  uncited  US  Army  technical 
reports  (Armstrong  et  al.,  1923;  Beebe,  1924; 
Marshall,  1917  a,b;  McElroy  et  al.,  1943).  Finally, 
previously  derived  human  toxicity  estimates 
for  chlorine  were  also  reviewed  and  collected 
(NRC,  2004;  Perry  &  Articola,  1980;  Rijnmond 
Public  Authority,  1982;  ten  Berge  &  van  Heemst, 
1983;  Withers  &  Lees,  1985a,b). 

[Some  of  the  LCT50S  in  Table  2  do  not  exactly 
correspond  to  what  was  originally  reported. 
Probit  analysis  (Finney,  1971)  with  modem  com¬ 
puter  software  was  not  always  available,  and 
many  older  studies  had  to  use  graphical  methods 
(Litchfield  &  Whcoxon,  1949).  Where  the  origi¬ 
nal  response  data  were  available,  the  originally 
reported  LCT50  was  updated  via  probit  analysis.] 

Only  studies  where  mammals  were  ex¬ 
posed  for  a  defined  duration  and  then  removed 
from  the  chamber  for  observation  were  used  in 
the  subsequent  statistical  analysis  (See  "DATA 
ANALYSIS  AND  RESULTS").  Studies  that  ex¬ 
posed  the  animals  continuously  until  death 
were  omitted.  [The  latter  studies  generally  pro¬ 
duced  higher  LCT50S  relative  to  fixed  duration 
and  observed  studies.  An  example  of  an  expo- 
sme  to  death  study  is  Weedon  et  al.  (1940),  which 
is  cited  by  many  reviews  of  chlorine  toxicity 
(MHAP,  1988;  Marman,  2005,  NRC,  2004)  and 
was  not  included  in  the  dataset  for  this  study] 

Many  sources  give  incomplete  information 
about  what  their  "lethal"  values  truly  signify 
{i.e.  is  it  1%  lethality  or  50%  lethality?),  so  only 
LCT50S  were  used  in  the  present  study.  If  an 
LCT50  could  not  be  estimated  from  any  reported 
quantal  response  data,  that  study's  "lethal" 
value  was  not  used.  Also,  no  LCT50  was  in¬ 
cluded  unless  the  exposure  duration  was  given. 

Every  effort  was  made  to  verify  that  each 
LCT50  was  an  independent  measurement.  For 
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Table  2.  Experimental  IH  LCT50S  and  Other  Parameter  Values. 


Study  (Authors) 

Report 

Date 

Raw 

Data 

Method 

Probit 

Slope 

Slope 

SE 

Time 

(min) 

LCT50 

(mg-min/ m®) 

LD50 

(mg) 

Mouse 

Silver,  McGrath 

1942 

Yes 

Probit 

7.32 

0.52 

10 

13150 

0.41 

Silver,  McGrath,  Ferguson 

1942 

Yes 

Probit 

7.83 

1.05 

10 

19430 

0.60 

Zwart,  Woutersen 

1988 

Yes 

Probit 

10.04 

2.01 

10 

30620 

0.95 

Alarie 

1980 

No 

10 

8752 

0.27 

Bunting 

1945 

No 

10 

15200 

0.47 

Lipton,  Rotariu 

1941 

No 

10 

18199 

0.56 

Bitron,  Aharonson 

1978 

Yes 

MLE 

10.63 

1.02 

11 

9245 

0.29 

Schlagbauer,  Henschler 

1967 

Yes 

Probit 

6.18 

1.23 

30 

10682 

0.33 

Zwart,  Woutersen 

1988 

Yes 

Probit 

6.82 

3.68 

30 

49300 

1.53 

Bitron,  Aharonson 

1978 

Yes 

MLE 

10.63 

1.02 

55 

27096 

0.84 

Back  et.  al.  &  Vernot  et.  al. 

1972  &  1977 

No 

60 

23822 

0.74 

Rat 

Zwart,  Woutersen 

1988 

Yes 

Probit 

15.84 

3.33 

5 

77850 

15.80 

Zwart,  Woutersen 

1988 

Yes 

Probit 

15.84 

3.33 

10 

61160 

12.42 

Zwart,  Woutersen 

1988 

Yes 

Probit 

15.84 

3.33 

30 

54810 

11.13 

Zwart,  Woutersen 

1988 

Yes 

Probit 

15.84 

3.33 

60 

83820 

17.02 

Vernot  et.  al. 

1977 

No 

60 

50947 

10.34 

Guinea  Pig 

Lehmann 

1882 

Yes 

MLE 

30 

290000 

52.78 

Rabbit 

Marshall 

1917(a) 

Yes 

MLE 

30 

21580 

31.36 

Barrow,  Smith 

1975 

Yes 

MLE 

30 

30479 

44.29 

Lehmann 

1882 

Yes 

MLE 

190 

300000 

435.90 

Cat 

Lehmann 

1882 

Yes 

MLE 

235 

390000 

329.55 

Dog 

Armstrong 

1923 

Yes 

Probit 

7.51 

2.83 

10 

58340 

246.78 

Marshall 

1917(a) 

Yes 

Probit 

3.12 

1.45 

30 

47010 

198.85 

Underhill 

1920 

Yes 

Probit 

3.45 

0.81 

30 

49380 

208.88 

Marshall 

1917(b) 

No 

30 

75000 

317.25 

Beebe 

1924 

No 

30 

90000 

380.70 

Sheep 

Batchinsky  et  al. 

2006 

No 

30 

24343 

338.37 

Goat 

McElroy  et  al. 

1943 

Yes 

MLE 

18.5 

51500 

494.40 

McElroy  et  al. 

1943 

Yes 

MLE 

8.5 

16500 

158.40 

instance,  both  Back  et  al.  (1972)  and  Vernot  et  al. 
(1977)  report  the  same  data  for  mouse  (attrib¬ 
uted  to  Back  et  al.  in  Table  2)  and  rat  (attributed 
to  Vernot  et  al.  in  Table  2). 

DATA  REDUCTION 

Prior  to  use,  several  calculation  procedures 
were  used  to  reduce  the  collected  data.  First, 
when  possible,  probit  analysis  was  used  to  esti¬ 
mate  both  kc  and  the  LCT50  (via  the  binary  logis¬ 


tic  regression  routine  in  MINITAB®).  If  a  study 
included  more  than  one  exposure  duration,  du¬ 
ration  was  treated  as  a  factor  (not  a  covariate)  in 
the  calculations  (see  "Inhalation  Dose-Response 
Statistics  and  the  Toxic  Load"),  thereby  produc¬ 
ing  a  best  estimate  for  the  true  slope  for  the 
whole  study  (which  was  the  case  with  Zwart 
&  Woutersen,  1988).  However,  in  many  cases, 
only  LCT50S  (and  no  other  data)  were  given. 

For  four  studies  (Barrow  &  Smith,  1975; 
Lehmarm,  1882;  Marshall,  1917a;  McElroy,  1943) 
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basically  only  limited  quantal  response  data 
were  recorded;  which  was  not  sufficient  for 
use  in  a  traditional  probit  analysis.  No  kc  or 
LCT50S  were  originally  reported.  For  these  sets, 
a  one-factor  MLE  analysis  (see  Fox,  1997  or 
Hulet  et  al.,  2006  for  examples)  with  a  probit  link 
function  was  used  to  estimate  the  LCT50. 

After  collection  or  estimation,  the  LCT50S 
were  converted  into  nominal  LD50S  using  Equa¬ 
tion  [4]  and  the  species  appropriate  M  and  Vm 
(expt.  avg.)  values  from  Table  1.  All  the  values 
shown  in  Table  2  were  used  in  the  final  model 
fit  of  Equation  [8].  The  LCT50S,  as  well  as  the 
previous  general  population  estimates  listed  in 
"Previous  Human  Lethality  Estimates  for  Chlo¬ 
rine  Inhalation  Toxicity",  are  shown  in  Eigure  2. 
The  LCT50  versus  Time  relationships  for  the 
Rijnmond  Report  (1982)  and  ten  Berge  &  van 
Heemst  (1983)  are  essentially  identical,  al¬ 
though  the  ten  Berge  &  van  Heemst  curve  is 
slightly  higher. 

DATA  ANALYSIS  AND  RESULTS 
Probit  Slope  Estimation 

Slope  Estimation  via  a  Weighted  Average.  Aweighted 
average  (Speigel,  1961)  was  taken  of  the  experi¬ 
mental  probit  slopes  listed  in  Table  2.  The  in¬ 
verses  of  the  variances  (the  squares  of  the 
standard  errors  (SE))  of  the  probit  slope  esti¬ 
mates  were  used  for  the  weights.  The  probit 
slopes  from  two  WWI  era  studies  (Marshall, 
1917a;  Underhill,  1920)  were  not  used  due  to 
concerns  about  the  quality  of  their  probit  slopes. 
It  was  found  that  the  weighted  average  equals 
7.96  with  95%  confidence  interval  of  7.21  to 
8.71.  It  is  assumed  that  the  above  weighted  aver¬ 
age  for  the  probit  slope  is  a  reasonable  estimate 
for  the  probit  slope  for  a  healthy  (i.e.  military) 
human  subpopulation. 

Slope  Estimation  from  Other  Sources.  Several  probit 
equations  for  human  lethality  have  been  devel¬ 
oped  or  estimated  within  the  past  30  years.  Five 
such  relations  are  listed  in  Table  3.  All  were  de¬ 
veloped  for  use  in  modeling  chlorine  toxicity 
within  the  general  population;  although  in  the 
case  of  Withers  &  Lees  (1985a,b),  no  probit  rela¬ 
tion  was  originally  reported  for  the  general  (or 
their  average)  population.  Instead,  they  reported 


estimated  concentrations  for  10,  50  and  90%  le¬ 
thality  for  two  separate  exposure  durations.  To 
derive  their  equation  [see  Table  3],  a  fit  was 
made  to  these  points.  Withers  &  Lees  also  devel¬ 
oped  probit  equations  for  the  healthy  and  vul¬ 
nerable  subpopulations. 

Four  of  the  five  equations  were  originally 
reported  using  probit  units,  natural  logarithms 
and  concentrations  expressed  in  ppm.  For 
ECBC  risk  assessment  and  toxicological  appli¬ 
cations,  normit  units,  base  10  logarithms  and 
concentrations  expressed  in  mg/ m^  are  typically 
used.  For  chlorine,  the  conversion  is  [mg  /  m^]  = 
[ppm]  X  (2.898).  Table  4  presents  the  probit  slopes 
(kc)  for  the  references  given  in  Table  3. 

Linear  Regression  Analysis 

Using  Equation  [8]  and  IRES  with  MINTTAB®, 
the  following  model  fit  was  obtained: 

log(LD5o)  =  (0.515)  +  (0.869)log(M) 

-f  (0.630)log(T) 

A  plot  of  Equation  [9]  (with  associated  95%  con¬ 
fidence  limits)  is  shown  in  Figure  3  with  the  in¬ 
dividual  experimental  LD50S  for  comparison. 

For  the  IRLS  process,  an  initial  unweighted 
model  fit  to  Equation  [8]  was  performed  using 
MINITAB®;  studentized  deleted  residuals  were 
calculated  for  the  set  of  observations.  From 
these  residuals,  a  set  of  weights  was  created  for 
the  next  iterahon.  Iterations  continued  unhl  a 
specific  tolerance  level  was  reached. 

Weights  (W)  were  calculated  as  follows 
(based  on  recommendations  of  Box  &  Draper, 
1987): 


convergence  for  each]  :  (W/  +  i,y  —  Wj-y) 
criteria  <1x10^ 

[10] 

where  the  counter  "i"  represents  ith  iteration, 
and  commences  at  2  since  the  first  regression 
run  is  unweighted,  “f  is  the  index  for  the  7*  ob¬ 
servation,  and  TRESi  j  is  the  studentized  deleted 
residual  for  the/^*’  observation  within  the  i*  iter¬ 
ation.  Convergence  was  defined  as  having  been 
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Figure  2.  Animal  Data  and  Previous  Human  Estimates  for  Chlorine  IH  Lethality  as  a  Function  of  Exposure  Duration. 
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Table  3.  Equations  for  Human  Chlorine  Inhalation  LCT50S. 

Reference  Year  Original  Published  Equation  *  Converted  Equation** 


Withers,  Lees 

1985 

Eisenberg, 

Lynch, 

Breeding 

1975 

Yp  (Probit)  =  -17.10  +  1.69 
Ln  (  C  T) 

Perry, 

Articola 

1980 

Yp  (Probit)  =  -36.45  +  3.13 
Ln  (  C  T) 

ten  Berge, 
van  Heemst 

1983 

Yp  (Probit)  =  -5.04  +  0.50 
Ln  (  C  T) 

Rijnmond 

Report 

(Harris, 

Moses) 

1982 

Yp  (Probit)  =  -11.40  +  0.82 
Ln  (  C  T) 

Ym  (Normit)  =  -11.697  +  1.675  LOG  (  C  T) 

Yn  (Normit)  =  -27.045  +  3.891  LOG  (  C  T) 

Yjv  (Normit)  =  -50.242  +  7.207  LOG  (  C  T) 

Yn  (Normit)  =  -11.503  +  1.151  LOG  (  C  T) 

Ym  (Normit)  =  -18.799  +  1.888  LOG  (  C  T) 


*  The  concentration  "C"  is  **  The  concentration  "C"  is  measured  in 

measured  in  parts  per  million  milligrams  per  cubic  meter 


reached  when  for  each  Wj  successive  iterations 
(i  and  i  +  1)  produced  estimates  that  agreed  to 
the  sixth  decimal  place. 

Chlorine  toxicity  correlates  (with  p  < 
0.001)  with  species  body  mass  (M).  On  an  abso¬ 
lute  nominal  dose  basis,  LD50  scales  to 
Thus,  for  larger  animals  a  greater  agent  dose — 
on  an  absolute  basis,  is  required  for  lethality. 
However,  on  a  per  kg  of  body  mass  basis,  larger 
mammals  are  more  sensitive  to  chlorine  than 
smaller  mammals  [LD50  (mg/kg)  scales  toM 

Toxic  Load  Exponent  and 
Time-Concentration  Relationship 

Estimation  from  Aliometric  Modeiing  of  Historicai  Data¬ 
base.  Based  on  Equation  [8],  n  equals  [1  /  (1  - 
fcj)].  From  the  LD50  fit  (Equation  [9]),  fcr  equals 
(0.630);  thus,  the  TEE  equals  2.70.  For  a  healthy 
human,  the  minute-volume  for  light  activity  is 
assumed  to  equal  0.015  m^/min,  and  the  mili¬ 
tary  standard  mass  is  assumed  to  be  70  kilo¬ 
grams  (US  Army  Chemical  School,  2005). 
Using  the  above  information  and  Equation  [7], 
the  following  relation  is  obtained  (shown  also 
in  Figure  4): 

C^^-^°)r  =  4.39X10^“  [11] 

The  standard  error  (SE  =  0.059)  for  the  fitted  co¬ 
efficient,  kt  from  Equation  [9]  can  be  used  to  pro¬ 
duce  approximated  95%  confidence  intervals 
for  n.  Random  values  were  generated  by  us¬ 


ing  a  normal  distribution  with  mean  equal  to 
the  fitted  k^  and  a  standard  deviation  equal  to 
the  SE  of  the  fitted  kf.  Using  the  relationship, 
[1  /  (1  -  kj)],  the  corresponding  random  n 
values  were  calculated.  The  95%  confidence  in¬ 
tervals  correspond  to  the  2.5  and  97.5  percen¬ 
tiles  of  the  random  n  values,  for  a  confidence 
interval  of  2.06  to  3.96  in  this  case.  This  low  pre¬ 
cision  is  possibly  due  to  the  relatively  narrow 
range  of  exposure  durations  represented  in  the 
historical  data,  with  the  bulk  of  them  ranging 
from  10  to  60  minutes.  The  presence  of  random 
error  between  the  various  studies  could  also  de¬ 
crease  the  precision  in  the  TEE  estimate.  Thus, 
another  approach  was  investigated  to  develop 
a  more  precise  estimate  for  the  TEE  (see  below). 

Estimation  from  Individual  Studies.  Only  two  studies 
listed  in  Table  2  involved  exposures  at  several 
durations:  Bitron  &  Aharonson  (1978)  and 
Zwart  &  Woutersen  (1988). 


Table  4.  General  Population  Probit  Slopes  [(Base  10) 
on  Concentration]. 


Reference 

Probit  Slope 
(with  respect 
to  Concentration) 

Withers,  Lees 

3.3 

Eisenberg  et  al. 

10.7 

Perry,  Articola 

19.0 

ten  Berge,  van  Heemst 

3.2 

Rijnmond  Report 

5.2 

Military  Operations  Research,  V15  N3  2010 


Page  69 


REVIEW  AND  ASSESSMENT  OF  CHLORINE  MAMMALIAN  LETHALITY  DATA  AND  THE 
DEVELOPMENT  OF  A  HUMAN  ESTIMATE 


10000 


(0 

c 

1 

o 


1000 


O) 

E 

_s  100 


Human  Estimate 
(Present  Work) 
95%  Conf.  Limits 

0 

Mouse 

Rat 

X 

B 

Guinea  Pig 

Rabbit 

Cat 

B 

Dog 

s 

Goat 

Sheep 

J _ I _ I  I  I  I 


Time  (minutes) 

Figure  3.  Allometric  Model  Fit  of  Nominal  IH  LD50  for  Healthy  Human  Adults  Exposed  to  Chlorine. 


Bitron  &  Aharonson  (1978)  exposed  476 
male  mice  at  two  constant  vapor  concentrations 
over  a  range  of  exposure  durations  [13  to  160 
minutes  for  170  ppm  (16  runs),  and  6  to  30  min¬ 
utes  for  290  ppm  (18  rums)].  The  dataset  is  well 
balanced,  with  a  good  distribution  of  percent  re¬ 
sponse  values  that  cover  at  least  5  to  95%  range 
for  each  concentration.  A  value  for  n  was  not 
originally  reported,  but  ten  Berge  &  van  Heemst 
(1983)  subsequently  calculated  a  value  of  3.5 
(with  95%  Cl  of  2.5  to  4.5)  from  this  dataset. 

Zwart  &  Woutersen  (1988)  exposed  male 
and  female  mice  and  rats  to  several  chlorine 
concentrations  at  four  exposure  durations. 
Two  hundred  twenty  rats  were  used  for  expo¬ 
sures  of  5, 10,  30  and  60  minutes;  100  mice  were 
used  for  exposures  of  10  and  30  minutes.  For 
each  exposure,  10  animals  (five  of  each  gender) 
were  exposed  in  a  common  chamber.  Unfortu¬ 
nately,  the  response  results  were  not  broken 
down  by  gender.  [Gender  differences  often  exist 


in  the  responses  of  rodents  to  toxicants — see 
Benton  et  al.,  2007].  Also,  this  dataset  is  not  as 
well  balanced,  particularly  for  the  rat,  as  that 
of  Bitron  &  Aharonson  (1978).  It  is  interesting 
to  note  that  Mannan  (2005),  in  his  extensive  re¬ 
view  of  the  existing  chlorine  toxicity  data,  did 
not  mention  the  results  of  Zwart  &  Woutersen. 
Considering  the  relative  merits  of  the  two  stud¬ 
ies,  the  data  of  Bitron  &  Aharonson  are  superior 
to  that  of  Zwart  &  Woutersen  for  estimating  n. 

For  this  study,  response  results  as  a  function 
of  chlorine  vapor  concentration  and  duration 
were  extracted  from  the  graphical  representa¬ 
tions  shown  in  Figures  9  and  11  of  Bitron  & 
Aharonson  (1978).  Using  the  binary  logistic  re¬ 
gression  routine  in  MINITAB®,  a  multifactor 
probit  analysis  was  then  performed  on  the  data. 
It  was  found  that  the  TLE  equals  3.18  with  95% 
Cl  of  2.90  to  3.46.  This  confidence  interval  is  con¬ 
siderably  narrower  than  was  originally  calcu¬ 
lated  by  ten  Berge  &  van  Heemst  (1983). 
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Figure  4.  Comparison  of  Human  (Military)  Estimate  for  Chlorine  IH  Lethality  from  Allometric  Fit  to  Previous 
Human  Estimates  and  Experimental  Animal  Data. 


Human  Lethality  Estimates 
for  Chlorine  Inhalation 

Three  parameters  are  needed  to  mathemat¬ 
ically  define  a  final  inhalation  toxicity  estimate 
(as  represented  in  Equation  [7]:  (1)  a  probit  slope 
to  describe  the  range  of  individual  sensitivities 
to  the  toxicant  (kj-L)',  (2)  the  LCT50  (or  toxic  load) 
at  some  reference  duration  {kg  when  Yjv  is  set 
equal  to  zero);  and  a  TEE  (n)  to  describe  how 
the  toxicity  changes  with  duration. 

In  the  following  two  sections,  the  derivation 
of  the  final  human  estimates  is  described  for 
both  a  military  subpopulation  and  the  general 
population.  These  estimates  are  shown  in  Eig- 
ure  5  in  comparison  to  previous  estimates  and 
the  existing  mammalian  database.  The  final 
military  subpopulation  estimate  is  slightly  dif¬ 
ferent  from  what  was  derived  via  allometric 
scaling  (Equation  [11])  for  reasons  explained  in 
greater  detail  below. 


Military  Subpopulation 

Probit  Slope.  Eor  the  probit  slope,  a  weighted  av¬ 
erage  from  the  available  data  (see  Table  2) 
yielded  an  estimate  of  8.0  (rounded  from  7.96), 
with  95%  Cl  of  7.2  to  8.7.  It  is  assumed  that  this 
estimate  is  appropriate  for  a  healthy  (military) 
human  subpopulation. 

Reference  Median  Lethal  Dosage.  Eor  the  median  le¬ 
thal  toxic  load,  there  is  some  leeway  for  a  recom¬ 
mended  value  due  to  the  large  amount  of  random 
error  present  among  the  various  experimental 
studies  (see  Table  2).  In  Eigure  4,  the  allometric 
fit  (Equation  [9])  is  shown  in  comparison  to  esti¬ 
mates  previously  derived  by  other  researchers. 
The  allometric  fit  of  the  estimates  derived  herein 
is  higher  than  that  for  previous  estimates,  with 
those  of  ten  Berge  &  van  Heemst  (1983)  and  the 
Rijnmond  Report  (1982)  being  closest  to  these 
revised  estimates.  Because  of  faulty  assumptions 


Military  Operations  Research,  V15  N3  2010 


Page  71 


(jUi/uiui-Bui)  “ion  , 


REVIEW  AND  ASSESSMENT  OF  CHLORINE  MAMMALIAN  LETHALITY  DATA  AND  THE 
DEVELOPMENT  OF  A  HUMAN  ESTIMATE 


1e+6 


1e+5 


1e+4 


1e+3 


10  100 

Time  (minutes) 

a.  LCT50  versus  Exposure  Duration  (T) 


10000 


1000 


100 


Time  (minutes) 

b.  LC50  versus  Exposure  Duration  (7) 

Figure  5.  Comparison  of  Final  Human  Estimates  Chlorine  IH  Lethality  with  Previous  Human  Estimates  and 
Experimental  Data. 


Page  72 


Military  Operations  Research,  V15  N3  2010 


REVIEW  AND  ASSESSMENT  OF  CHLORINE  MAMMALIAN  LETHALITY  DATA  AND  THE 

DEVELOPMENT  OF  A  HUMAN  ESTIMATE 


used  in  their  derivations,  the  estimates  of 
Eisenberg  et  al.  (1975)  and  Perry  &  Articola 
(1980)  can  be  dropped  from  consideration. 

There  are  two  reasons  for  the  difference  be¬ 
tween  Equation  [9]  and  ten  Berge  &  van  Heemst 
and  the  Rijnmond  Report.  First,  Equahon  [9]  is 
meant  for  use  with  a  healthy  (military)  subpopu¬ 
lation,  while  the  other  two  are  for  the  general 
population.  [This  issue  is  further  addressed  be¬ 
low.]  Second,  Equation  [9]  is  based  on  the  fact 
that  both  the  LD50  (absolute)  and  LC50  increase 
as  a  function  of  increasing  species  body  mass 
(as  demonstrated  by  statistical  analysis  in  the 
present  study),  but  both  ten  Berge  &  van  Heemst 
and  the  Rijnmond  Report  assumed  that  a  human 
estimate  can  be  modeled  on  an  average  for  all 
mammals.  So,  as  a  starting  point  for  a  human 
general  population  estimate,  the  human  healthy 
(military)  LCT50  for  a  two  minute  exposure  was 
taken  from  Equation  [9],  or  13,500  mg-min/m^. 

Toxic  Load  Exponent.  For  deriving  the  TLE,  past 
approaches  have  been  somewhat  subjective, 
ten  Berge  &  van  Heemst  (1983)  took  an  average 
of  the  lethality  study  of  Bitron  &  Aharonson 
(1978)  (n  =  3.5)  and  the  olfactory  irritation  study 
of  Anglen  (1981)  (n  =  1.9)  to  arrive  at  2.75.  This 
may  be  a  case  of  obtaining  a  correct  answer  via 
incorrect  methods:  it  is  unlikely  that  the  toxic 
mechanisms  and  their  corresponding  n  values 
are  the  same  for  lethality  and  nuisance  effects. 
Thus,  an  average  of  the  two  is  not  the  proper  ap¬ 
proach  for  obtaining  a  lethal  TLE. 

Harris  &  Moses  (1983)  concluded  that  the 
Coast  Guard  reports  (Eisenberg  et  al.,  1975; 
Perry  &  Articola,  1980)  have  the  TLE  correct  at 
2.75.  It  is  not  readily  evident  how  the  Coast 
Guard  arrived  at  this  estimate,  but  it  appears 
that  instead  of  basing  it  on  the  available  mam¬ 
malian  data,  they  performed  an  empirical  fit  of 
previous  human  estimates. 

The  MHAP  (1988)  also  concluded  that  2.75 
was  a  reasonable  estimate  for  the  TLE.  They  rea¬ 
soned  that  the  evidence  for  this  value  is  much 
stronger  than  for  the  lower  values  derived  by 
others  (such  as  Withers  &  Lees,  1985a,b).  They 
believed  that  the  value  of  2.75  is  advantageous 
because  it  results  in  conservative  assessments 
for  short  exposures  (less  than  ten  minutes),  for 
which  there  is  considerable  uncertainty  in  the 
dose-response  relationship. 


Withers  &  Lees  (1985a)  based  their  estimate 
(n  =  2)  on  a  review  of  Weedon  et  al.  (1940)  and 
Bitron  &  Aharonson  (1978).  Others  (Franks 
et  al.,  1996;  Turner  &  Fairhurst,  1990)  have 
used  the  same  value-hut  for  different  reasons. 
Withers  &  Lees  felt  that  these  were  the  only 
available  studies  from  which  to  estimate  n  for 
chlorine.  A  major  difference  between  Weedon 
et  al.  and  Bitron  &  Aharonson  is  the  number  of 
animals  used  in  each  study,  24  and  476,  respec¬ 
tively:  In  essence.  Withers  &  Lees  took  an  un¬ 
weighted  average  of  the  n  values  from  these 
studies.  All  other  factors  being  equal,  more 
weight  should  have  been  given  to  the  Bitron  & 
Aharonson  value,  since  it  is  based  upon  a  larger 
number  of  animals.  If  a  weighted  average  is 
taken  (using  the  number  of  animals  for  the 
weights),  an  n  of  2.94  is  obtained. 

In  summary,  this  study  is  in  agreement  with 
both  ten  Berge  &  van  Heemst  (1983)  and  Withers 
&  Lees  (1985a)  on  the  importance  of  Bitron 
&  Aharonson  (1978).  The  results  of  Bitron  & 
Aharonson  support  a  higher  value  (in  the  range 
of  3  to  3.5)  for  the  TLE  than  has  been  used  pre¬ 
viously.  Of  the  available  experimental  studies, 
Bitron  &  Aharonson  have  the  best  set  of  lethality 
data  involving  multiple  durations,  but  this  is 
only  one  study  which  has  not  yet  been  repro¬ 
duced.  As  a  result,  cauhon  is  required  with  re¬ 
spect  to  TLE  recommendation.  Setting  the  TLE 
value  too  high  will  underestimate  the  toxic  ef¬ 
fects  for  longer  durations.  The  range  of  reason¬ 
able  values  to  choose  from  can  be  defined  on 
the  low  end  by  the  allometric  fit  from  this  study 
(n  =  2.70  from  Equahon  [9])  and  on  the  upper 
end  by  this  study's  estimate  of  n  from  the  data 
of  Bitron  &  Aharonson  (n  =  3.18).  Thus,  for  a  hu¬ 
man  estimate,  a  TLE  of  2.75  was  chosen  (as  a  con¬ 
servative  estimate)  for  generating  a  final  human 
(military)  toxicity- time  relationship.  This  favors 
the  bulk  of  the  historical  database  (Table  2)  and 
previous  estimates  for  the  exponent  value  made 
by  previous  researchers  (Eisenberg  et  al.,  1975; 
Harris  &  Moses,  1983;  MHAP,  1988;  Perry  & 
Arhcola,  1980;  Rijnmond  Public  Authority,  1982; 
ten  Berge  &  van  Heemst,  1983).  However,  this 
TLE  should  be  re-evaluated  in  the  light  of  any  ap¬ 
plicable  future  mammalian  studies. 

Human  (Military)  Lethality  Relationship  from  Allometric 
Fit.  From  the  previous  discussion,  the  final  human 
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toxicity  estimate  for  healthy  (military)  individ¬ 
uals,  based  upon  allometric  scaling,  is  defined 
by  Equation  [12]  (see  Figure  5a): 

Yn  =  (  -31.52)  +  (2.91)log(C^-^®T) 

or  [12] 

L(C^-^®T)5o  =  6.79  X  10^° 
where  C  is  in  mg/ m^  and  T  is  in  minutes. 

General  Population.  Toxicity  estimates  for  the  gen¬ 
eral  population  were  generated  from  the  esti¬ 
mate  for  the  military  subpopulation  (Equation 
[12])  via  the  method  of  Crosier  (2007)  (as  previ¬ 
ously  discussed).  Using  this  method,  the  follow¬ 
ing  values  were  obtained: 

Probit  slope  (general  population)  (kc)  =  6.0 

LCT50  for  two  minutes  (general  popula¬ 
tion)  =  9,500  mg-min/m^ 

For  the  TEE,  mathematical  considerations  re¬ 
quire  that  both  a  subpopulation  and  the  general 
population  must  have  identical  exponent  values 
{ex.  differing  values  for  n  can  produce  the  par¬ 
adoxical  situation  where  the  LCT50  for  the 
general  population  could  exceed  that  of  the  mil¬ 
itary  subpopulation).  Putting  these  three  pa¬ 
rameters  together  produces  the  probit  relation 
for  the  general  population: 

Yn  =  (  -22.698)  +  (2.18)log(C^-^®T) 

or  [13] 

L(C^-^®T)5o  =  2.58  X  10“ 
where  C  is  in  mg/ m^  and  T  is  in  minutes. 

HPAC  Modeling  Results 

A  series  of  modeling  runs  was  conducted, 
using  the  Defense  Threat  Reduction  Agency's 
(DTRA)  Hazard  Prediction  and  Assessment  Ca¬ 
pability  (HPAC,  2003).  For  a  typical  large-scale 
accidental  release  of  chlorine,  vapor  concentra¬ 
tion  versus  time  profiles  were  generated  for  sev¬ 
eral  locations  downwind.  Then  several  human 
(general  population)  toxicity  estimates  were 
used  to  estimate  the  probability  of  an  exposed 
individual  dying  should  the  person  be  exposed 
to  these  concentration-time  profiles. 

HPAC  Model  Inputs  and  Outputs.  The  scenario  mod¬ 
eled  was  for  a  catastrophic  release  of  50  tons 


(—45,500  kg)  of  chlorine  liquid.  This  is  based 
on  the  largest  release  reported  by  Mannan 
(2005)  in  his  review  of  chlorine  industrial/ trans¬ 
portation  accidents. 

Model  inputs  were  kept  fairly  simple.  Flat  ter¬ 
rain  was  assumed.  Meteorological  condihons 
were  kept  at  a  steady-state,  fixed  wind  direction 
and  speed  during  the  simulations.  The  HPAC 
model  estimates  the  amormt  of  chlorine  that  ini¬ 
tially  pools,  but  aU  the  liquid  evaporated  very 
quickly.  Since  meteorological  conditions  can  pro¬ 
duce  significant  differences  in  hazard  prediction 
results,  three  different  weather  conditions  were 
modeled: 

•  Low  wind  (1  m/sec),  clear  sky,  nighttime 
(2am),  "stable"  condition  (Pasquill  Stability 
(PST)  Class  F— Pasquill,  1961) 

•  Moderate  wind  (5  m/  sec),  cloudy  day,  "neu¬ 
tral"  condition  (PST  Class  D) 

•  Low  wind  (2  m/sec),  clear  sky,  daytime 
(2pm),  "unstable"  condition  (PST  Class  B). 

The  HPAC  model  provided  numeric  output  of 
concentration  time  history  for  each  model 
run  at  downwind  distances  of  0.2,  0.5,  1.0,  1.5 
and  2.0  km.  The  concentration  time  histories 
were  then  numerically  integrated  using  several 
different  human  toxicity  estimates  from  Table  3 
and  Equation  [13]  to  calculate  the  probability  of 
lethality  at  each  downwind  sampler  location. 

Samples  of  some  of  the  concentration-time 
and  TL-time  outputs  as  a  function  of  downwind 
distance  are  shown  in  Figures  6a  and  6b  for  the 
PST  Class  D  HPAC  runs.  The  peak  concentration 
values  were  observed  rather  quickly,  from  es- 
senhally  instantaneous  at  0.2  km  to  about  7  min¬ 
utes  for  2  kms  downwind.  The  peaks  roughly 
correspond  to  the  point  where  the  accumulated 
TL  had  reached  about  90  to  95%  of  its  final  value 
for  a  particular  downwind  location.  If  the  accu¬ 
mulated  TL  had  not  reached  lethal  levels  within 
a  few  minutes  after  the  observed  peak  concen¬ 
tration,  it  was  very  unlikely  that  it  will  do  so. 
Of  the  five  downwind  locations,  only  0.2  and 
0.5  km  had  an  accumulated  TL  that  was  in  ex¬ 
cess  of  the  LTL05  from  either  this  study  (Equa¬ 
tion  [12])  or  the  Rijnmond  Report  (see  Table  3). 

Sensitivity  of  HPAC  Model  Predictions  to  Changes  in 
Toxicity  Estimate  Parameters.  Not  surprisingly,  as 
demonstrated  in  Figure  7  using  the  HPAC 
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a.  Concentration  (2  min  TWA)  versus  Time  Profile 


Time  (minutes) 

b.  Toxic  Load  versus  Time  Profile 

Figure  6.  Output  Profiles  from  HPAC  as  a  Function  of  Downwind  Distance  for  Chlorine  (Pasquill  Stability 
Category  D). 
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— A — 

LTLjj  =  2.582  x  10’°,  =  6.0  &  n  =  2.75 

-  -V  - 

LTL5,  =  1.267  X  10',  kj  =  6.0  &  n  =  1.85 

— -D— 

LTLsj  =  2.582  x  10’°,  k,.  =  3.0  8.  n  =  2.75 

— 0- 

LTLj,  =  1.267  X  10',  kc  =  3.0  &  n  =  1.85 

Downwind  Distance  (km) 

Figure  7.  Probability  of  Lethality  versus  Downwind  Distance  as  a  Function  of  Toxic  Load  Expression  (Pasquill 
Stability  Category  F). 


concentration-time  profile  generated  under 
night  time  conditions  (PST  Class  F),  the  relation¬ 
ship  of  the  probability  of  lethality  as  a  function 
of  downwind  distance  is  dependent  on  the  tox¬ 
icity  estimate  used.  Four  curves  were  generated 
using  a  two-by-two  matrix  of  kc  (3  and  6)  and  n 
(1.85  and  2.75)  values.  All  four  curves  share  an 
identical  LC50  of  4,750  mg/m^  at  two  minutes. 
Equation  [13]  (for  the  general  population)  is 
represented  by  the  curve  labeled  =  6  and 
n  =  2.75. 

An  increase  in  the  TLE  value  shifts  the  line 
to  the  left,  with  the  new  curve  still  parallel  to 
the  original  line  [when  plotted  on  a  probability 
versus  log(distance)  scale].  An  increase  in  kc  ro¬ 
tates  the  new  line  clockwise  relative  to  the  orig¬ 
inal  line.  Changing  both  kc  and  n  both  shifts  and 
rotates  the  new  line  compared  to  the  original. 
Since  a  probit  slope  change  rotates  the  line,  the 
differences  in  downwind  distances  between 


the  old  and  new  lines  is  a  function  of  the  lethal¬ 
ity  probability,  with  the  difference  increasing  as 
one  moves  closer  or  further  from  the  point  of  in¬ 
tersection  between  the  two  lines.  These  trends 
from  Figure  7  are  better  highlighted  in  Figure  8. 

Comparison  of  Human  Toxicity  Estimates  in  Predicting 
Lethaiity  Probabilities  from  HPAC  Model  Predictions. 
In  Figure  9,  the  probability-distance  profile 
based  upon  Equation  [13]  is  compared  to  pro¬ 
files  generated  using  the  human  toxicity  esti¬ 
mates  in  Table  3  for  the  HPAC  concentration¬ 
time  profile  generated  for  night  time  conditions 
(PST  Class  F).  For  probabilities  below  10%,  the 
downwind  distance  predicted  by  Equation  [13] 
is  less  than  what  is  predicted  by  Withers  & 
Lees  (1985b),  Rijnmond  Report  (1982)/Harris  & 
Moses  (1983)  and  ten  Berge  &  van  Heemst 
(1983).  Distances  predicted  by  Haber's  Rule  (us¬ 
ing  a  LCT50  of  9,300  mg-min/m^  and  a  probif 
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Downwind  Distance  (km)  Downwind  Distance  (km) 


Change  in  slope  rotates  plot 


Change  in  TLE  shifts  plot 


Figure  8.  Trends  in  Lethality  Probability  versus  Downwind  Distance  as  a  Function  of  Probit  Slope  (Concentra¬ 
tion)  and  Toxic  Load  Exponent  (TLE). 


slope  of  6)  were  also  calculated.  Considering  the 
historical  record  (see  HISTORICAL  RECORD 
OF  CHLORINE  RELEASES),  the  downwind 
distances  predicted  by  Haber's  rule  are  very  un¬ 
realistic  {i.e.  a  20%  chance  of  lethality  at  a  down¬ 
wind  distance  of  2  km). 

Of  the  other  three  studies,  the  Rijnmond  Re¬ 
port  curve  is  the  closest  in  location  and  orienta¬ 
tion  to  that  of  Equation  [13].  The  probit  slope 
and  TLE  values  from  the  Rijnmond  Report 
(k-ri  =  1.89  and  n  =  2.75)  are  very  similar  to  that 
of  Equation  [13]  (kri  =  2.18  and  n  =  2.75).  The 
difference  between  the  two  curves  is  mainly 
due  to  differences  in  the  LCT50  used  in  the  two 
estimates.  Both  Withers  &  Lees  and  ten  Berge  & 
van  Heemst  estimated  smaller  probit  slopes 
than  Equation  [13]  and  the  Rijnmond  Report, 
and  their  prediction  curves  are,  thus,  rotated 
counter-clockwise  relative  to  Equation  [13]  and 
the  Rijnmond  Report.  In  Figure  10,  the  same 
toxicity  relationships  used  for  Figure  9  are  used 
for  the  HPAC  concentration-time  profile  gener¬ 
ated  for  day  time  conditions  (PST  Class  B). 

In  Figure  11,  the  previous  plots  for  Equation 
[13],  ten  Berge  &  van  Heemst  and  the  Rijnmond 
Report  are  shown  for  both  day  and  night  time 
conditions  (plots  for  overcast  conditions- -PST 
Class  D)  are  very  similar  to  those  of  night  time). 
Any  proposed  toxicity  estimate/ATD  model 
combination  should  produce  downwind  pre¬ 
diction  curves  consistent  with  the  historical  re¬ 
cord  (maximum  fatality  downwind  distance 
of  400  meters).  Under  daytime  conditions,  the 
plots  for  Equation  [13],  ten  Berge  &  van  Heemst, 
and  the  Rijnmond  Report  are  consistent  with  the 


historical  record.  However,  for  nighttime  condi¬ 
tions,  all  three  toxicity  estimates  still  produce 
significant  lethality  probabilities  at  400  meters 
downwind  (from  40  to  75%).  From  this  point 
downwind,  probabilities  from  both  Equation 
[13]  and  the  Rijnmond  Report  fall  off  quickly 
as  a  function  of  distance  relative  to  that  pre¬ 
dicted  by  ten  Berge  &  van  Heemst.  It  is  possible 
that  deaths  have  not  been  seen  at  500  meters 
(despite  the  high  predicted  probabilities)  due 
to  other  factors  {i.e.  more  time  available  for  es- 
cape/sheltering  in  place  at  the  longer  down¬ 
wind  distances,  night  time  accidents  may  not 
be  that  frequent  relative  to  day  time,  not  all  re¬ 
leases  are  50  tons,  etc.).  At  1  km,  there  is  still 
a  3%  probability  of  lethality  with  ten  Berge  & 
van  Heemst,  while  the  other  two  are  both  below 
0.2%.  This  greater  distance,  on  the  part  of  ten 
Berge  &  van  Heemst,  is  due  primarily  to  its  shal¬ 
lower  probit  slope  {kc  =  3.2  versus  5.19  and  6  for 
the  other  two  toxicity  estimates). 


DISCUSSION/CONCLUSIONS 

The  development  of  a  human  toxicity  esti¬ 
mate  for  chlorine  has  been  hampered  by  the  lack 
of  appropriate  and  accurate  toxicity  data.  As  a 
result,  various  investigators  have  had  to  resort 
to  various  degrees  of  subjective  "hand-waving" 
to  produce  appropriate  estimates.  Some  (Harris 
&  Moses,  1983;  Richmond  Public  Authority, 
1982;  ten  Berge  &  van  Heemst,  1983;  Withers 
&  Lees,  1985a,b)  have  been  more  successful  than 
others  (Eisenberg  et  al.,  1975;  Perry  &  Articola, 
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Figure  9.  Probability  of  Lethality  versus  Downwind  Distance  for  Several  Human  Lethality  Estimates  for 
Chlorine  (Pasquill  Stability  Category  F). 


1980).  The  present  study  improved  upon  the 
previous  work  by  applying  novel  approaches 
{i.e.  allometric  modeling,  weighted  average  of 
experimental  probit  slopes,  statistical  analysis 
of  original  response  data,  mathematical  model¬ 
ing  of  relationship  between  subpopulation  and 
the  general  population,  etc.)  to  the  problem 
and  introducing  previously  unknown/under 
analyzed  lethality  data  (Armstrong  et  al.,  1923; 
Barrow  &  Smith,  1975;  Beebe,  1924;  Lehmann, 
1882;  Marshall,  1917a,b;  McElroy  et  al,  1943). 
Because  of  significant  shortcomings  in  the  chlo¬ 
rine  lethality  data,  one  should  not  be  dogmatic 
about  the  human  estimate.  Thus,  the  human 
toxicity  estimates  from  the  present  study  (for 
both  military  and  general  population)  are  of¬ 
fered  as  qualitative  (rather  than  quantitative) 
improvements  over  estimates  from  previous  in- 
vesfigators  (Harris  &  Moses,  1983;  Richmond 
Public  Aufhorify,  1982;  ten  Berge  &  van  Heemst, 
1983;  Withers  &  Lees,  1985a,b).  The  following  is 


a  brief  review  of  the  salient  issues  addressed  in 
the  present  study  and  how  they  differentiate  the 
recommendations  of  the  present  study  from 
those  of  the  previous  efforts. 


Variability  in  Human  Response 

It  has  been  recognized  by  many  that  the 
general  human  population  is  likely  to  have 
more  inhalation  response  variance  than  healthy 
laboratory  animals  (Crosier,  2003,  2007;  Crosier 
&  Sommerville,  2002;  Mannan,  2005;  MHAP, 
1988;  Richmond  Public  Authority,  1982;  Salem, 
1987;  Schubach,  1997;  Sommerville  et  al,  2005; 
ten  Berge  &  van  Heemst,  1983;  Withers  &  Lees, 
1985a,b).  The  problem  has  been  how  to  best  ex¬ 
trapolate  probit  slopes  from  a  laboratory  animal 
to  a  general  human  basis. 

In  the  present  study,  more  statistical  rigor 
was  employed.  This  was  done  by  using  probit 
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Figure  10.  Probability  of  Lethality  versus  Downwind  Distance  for  Several  Human  Lethality  Estimates  for 
Chlorine  (Pasquill  Stability  Category  B). 


analysis  to  generate  probit  slopes  and  associ¬ 
ated  standard  errors  from  the  original  experi¬ 
mental  data.  Then,  a  weighted  average  of  fhe 
probif  slopes  was  calculafed,  and  this  was  as¬ 
sumed  to  be  the  best  estimate  for  variabilify  in 
the  response  among  study  animals.  It  was  rea¬ 
sonably  assumed  that  healthy  humans  were 
equivalent  in  their  variability  to  laboratory  ani¬ 
mals.  An  estimate  for  fhe  general  population 
was  generated  by  using  a  mathematical  model 
of  fhe  relafionship  befween  a  subpopulation 
and  the  whole  population  developed  by  Crosier 
(2007).  Of  the  previous  estimates  for  fhe  probit 
slope  of  the  general  population,  the  estimate  of 
the  present  study  {kc  =  6)  is  in  closest  agreement 
to  that  of  fhe  Rijnmond  Reporf  {kc  =  5.2).  The 
other  important  studies  (ten  Berge  &  van 
Heemst  and  Withers  &  Lees)  arrived  at  lower 
probit  slopes  (indicating  more  response  vari¬ 
ability). 


Time  Dependence  of  Chlorine 
Inhalation  Toxicity 

Previous  researchers  (ten  Berge  &  van 
Heemst,  1983;  Withers  &  Lees,  1985a,b)  have 
recognized  the  futility  of  calculating  an  estimate 
for  fhe  TLE  based  on  the  existing  lethality  data¬ 
base.  Instead,  the  individual  findings  of  four 
experimental  toxicity  studies  (Anglen,  1981; 
Bitron  &  Aharonson,  1978;  Weedon  et  al.,  1940; 
Zwart  &  Woutersen,  1988)  were  cited  as  the  ba¬ 
sis  for  a  human  lethalify  exponent  value  by  sev¬ 
eral  studies  (ten  Berge  &  van  Heemst,  1983; 
Withers  &  Lees,  1985a,b;  Zwart  &  Woutersen, 
1988).  Although  their  work  has  been  misapplied 
in  the  past,  the  present  study  has  concluded  that 
only  Bitron  &  Aharonson  is  suitable  for  use  in 
developing  a  foxic  load  exponenf  for  lefhalify 
(see  Human  Lethality  Estimates  for  Chlorine 
Inhalation).  The  proper  use  of  Bifron  & 
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Downwind  Distance  (km) 

Figure  11.  Probability  of  Lethality  versus  Downwind  Distance  for  Several  Human  Lethality  Estimates  for 
Chlorine  (Pasquill  Sfability  Categories  B  and  F). 


Aharonson  leads  to  a  TLE  ranging  from  3  to  3.4. 
However,  allometric  modeling  and  consider¬ 
ation  of  the  database  as  a  whole  support  the 
slightly  lower  exponent  of  about  2.75.  The 
"truth"  probably  lies  somewhere  between  the 
two  extremes.  The  lower  exponent  (2.75)  is  rec¬ 
ommended  for  two  reasons:  (1)  there  is  no  collab- 
orahng  study  confirming  Bitron  &  Aharonson; 
and  (2)  it  is  desirable  to  avoid  underestimating 
toxicity  at  longer  durations  by  using  a  TLE  that 
is  too  large. 


Median  Lethal  Dosage  (LCT50) 

With  respect  to  LCT50,  the  basic  philosophy 
for  extrapolating  from  animals  to  humans  has 
been  to  assume  that  humans  are  equivalent  to 
the  average  mammal  (Harris  &  Moses,  1983; 
Rijnmond  Public  Authority,  1982;  ten  Berge  & 
van  Heemst,  1983;  Withers  &  Lees,  1985a,b). 


However,  although  the  exact  reason  has  been 
debated  (MHAP,  1988;  Withers  &  Lees,  1985a), 
it  has  been  recognized  that  the  chlorine  inhala¬ 
tion  LCT50  does  increase  as  body  mass  in¬ 
creases.  An  allometric  fit  of  the  data  in  the 
present  study  found  that  the  LCT50  was  slightly 
dependent  on  the  body  mass  (with  statistical 
significance),  with  both  the  nominal  LD50  (abso¬ 
lute)  and  LCT50  increasing  as  the  body  mass  in¬ 
creases.  However,  the  large  amount  of  random 
error  present  in  the  mammalian  database  (Table 
2)  decreases  the  confidence  in  allometric-based 
extrapolation  from  this  database  to  arrive  at 
a  human  estimate.  Thus,  the  "true"  answer 
probably  lies  in-between  the  two  philosophies. 

Final  Human  Lethality  Estimates 

Of  the  five  previous  human  estimates 
reviewed  (see  Table  3),  the  Rijnmond  Report 
(1982)  /Harris  &  Moses  (1983)  comes  the  closest 
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to  the  general  population  estimate  of  the  present 
work  (Equation  [13])  in  terms  of  downwind 
hazard  predictions  (see  HPAC  Modeling  Re¬ 
sults).  Both  have  identical  TLEs  (2.75)  and  sim¬ 
ilar  probit  slopes  (kc)  (6  versus  5.2). 

The  HPAC  modeling  runs  (see  HPAC  Mod¬ 
eling  Results)  confirms  the  often  observed  dis¬ 
connect  between  the  historical  record  (see 
HISTORICAL  RECORD  OF  CHLORINE  RE¬ 
LEASES)  and  the  overly  pessimistic  predictions 
from  ATD  models  (Mannan,  2005).  There  are 
three  possible  explanations  for  this  conflict.  Ei¬ 
ther,  the  ATD  models  are  over-predicting  the 
downwind  transport  of  chlorine;  the  human  le¬ 
thality  estimates  are  over-predicting  the  toxicity; 
or  a  combination  of  both  of  these  factors.  There  is 
a  temptation  to  be  biased  towards  using  larger 
LCT50  eshmates  to  achieve  more  reasonable  ca¬ 
sualty  predictions  from  the  current  ATD  models. 
However,  future  improvements  in  the  ATD 
models  (and/or  the  knowledge  of  the  atmo¬ 
spheric  chemistry  of  chlorine)  could  leave  such 
a  biased  lethality  eshmate  “out  of  posihon." 

Eor  the  three  important  toxicity  estimate  pa¬ 
rameters  (probit  slope  (kc),  TLE  (n),  and  refer¬ 
ence  LCT50),  the  general  population  estimate 
from  Equation  [13]  is  in  better  agreement  with 
the  Rijnmond  Report  than  with  the  other  previ¬ 
ous  estimates,  particularly  for  kc  and  n.  The  con¬ 
stant  difference  in  LCT50S  is  due  to  differences 
in  extrapolation  philosophies  (average  mammal 
(Rijnmond)  versus  allometric  scaling  (present 


study)).  The  LCT50  estimate  of  Equation  [13]  is 
about  a  factor  of  1.5  higher  than  that  from  the 
Rijnmond  Report,  which  is  not  much  above 
the  random  noise  in  the  data.  The  true  three  pa¬ 
rameter  fit  for  general  population  toxicity  is 
probably  bracketed  by  Equation  [13]  and  the 
Rijnmond  Report. 

In  summary,  the  following  upper  and  lower 
general  population  estimates  for  chlorine  lethal¬ 
ity  are  presented  for  consideration.  The  upper 
estimate  is  defined  by  Equation  [13],  and  the 
lower  estimate  is  defined  by  a  modified  form 
of  the  Rijnmond  Report  expression  (see  Table 
4).  The  Rijnmond  Report  expression  was  modi¬ 
fied  by  changing  the  original  probit  slope  (kc) 
from  its  value  of  5.2  equal  to  that  of  the  slope 
in  Equation  [13]  (kc  =  6.0),  and  then  a  two  min¬ 
ute  reference  dosage  of  6,500  mg-min/ m^  (cal¬ 
culated  from  the  Rijnmond  estimate)  was  used 
as  the  basis  for  a  toxic  load  expression  (with 
no  change  in  the  value  of  n).  In  comparison, 
the  two  minute  general  population  LCT50  from 
Equation  [13]  equals  9,500  mg-min/m^.  These 
estimates  are  shown  below  in  Equations  [13] 
and  [14]  and  in  tabular  form  in  Table  5. 


upper 

estimate 


Yn  =  (  -22.698)  +  (2.18)log(d-^^r) 


(adjusted 

allometric 


fit) 


or 

L(c2.75t)^^  =  2.58  X  10“ 


[13] 


Table  5.  General  and  Military  Populations  Dosage  and  Concentration  Estimates 


Modified  Rijnmond 

General  Population  General  Population  Military  Population 


Equation  [14]  Equation  [13]  Equation  [12] 


LCT50 
Time  (mg-min  / 

LC50 

m^)  (mg  /  m 

LC50  LCT50 

^)  (ppm)  (mg-min  / 

LC50 

m^)  (mg  /  m‘ 

LC50  LCT50 

’)  (ppm)  (mg-min  / 

LC50 

m^)  (mg  /  m 

LC50 

3)  (ppm) 

2 

6500 

3250 

1120 

9500 

4750 

1640 

13500 

6750 

2320 

10 

18100 

1810 

630 

26400 

2640 

910 

37600 

3760 

1300 

30 

36400 

1210 

420 

53200 

1780 

610 

75700 

2520 

870 

60 

56600 

940 

330 

82700 

1380 

480 

118000 

1960 

680 

120 

88000 

730 

250 

129000 

1080 

370 

183000 

1520 

530 

240 

137000 

570 

200 

200000 

830 

290 

284000 

1180 

410 

360 

177000 

490 

170 

259000 

720 

250 

368000 

1020 

359 

480 

213000 

440 

150 

311000 

650 

220 

442000 

920 

320 

960 

330000 

340 

120 

483000 

500 

170 

687000 

720 

250 
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lower  Yn  =  (  -21.710)  +  (2.18)log(d-^®r) 
estimate 
(modified 

Rijnmond)  L(C^'^^T)5q  =  9.09  X  10^ 

[14] 

To  evaluate  the  relative  merits  of  Equations  [13] 
and  [14],  it  is  recommended  that  a  more  exten¬ 
sive  sensitivity  study  with  ATD  models  be  con¬ 
ducted  (beyond  the  scope  of  the  present  work). 
If  ATD  models  continue  to  overstate  the  hazard 
(even  when  using  the  higher  estimate  of  Equa¬ 
tion  [13]),  then  it  is  very  likely  that  the  problem 
lies  with  the  algorithms  of  the  ATD  models 
(and/or  inadequate  knowledge  of  cloud  trans¬ 
port  or  chlorine  atmospheric  chemistry),  rather 
than  an  inadequate  lethality  estimate.  A  general 
population  median  lethality  estimate  higher 
than  Equation  [13]  cannot  be  justified  using 
the  current  lethality  database  (Table  2). 

RECOMMENDATIONS 

To  better  refine  the  human  lethality  esti¬ 
mates  for  chlorine  toxicity,  future  toxicity  stud¬ 
ies  should  address  the  following  issues:  Eirsf, 
there  is  currently  only  one  experimental  study 
(Bitron  &  Aharonson,  1978)  that  provides 
a  defensible  estimate  for  the  lethality  TEE  (n). 
The  reproducibility  of  fhis  result  should  be  con¬ 
firmed  (over  a  range  of  durations  from  several 
minutes  to  several  hours),  preferably  in  multi¬ 
ple  mammalian  species.  Second,  the  usefulness 
of  allomefric  scaling  can  be  strengthened  by 
choosing  larger  species  {i.e.  rabbit,  swine  or 
goat)  for  the  toxicity-time-dependency  studies. 
Lastly,  comprehensive  sensitivity  studies  using 
ATD  models  should  be  conducted  to  evaluate 
against  the  historical  record  casualty  predic¬ 
tions  from  all  toxicity  estimate  candidates. 
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HPAC  Hazard  Prediction  and  Assessment 
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IH  Inhalation 

IRES  Iteratively  reweighted  least  squares 

LCT50  Lethal  concentration-time  to  50%  of 

exposed  individuals 

LD50  Lethal  dosage  to  50%  of  exposed  indi¬ 
viduals 

LTL50  Lethal  toxic  load  50%  of  exposed  indi¬ 
viduals 

MHAP  Major  Hazards  Assessment  Panel 
MLE  Maximum  Likelihood  Estimation 
PST  Pasquill  Stability 

SE  Standard  error 

TL  Toxic  load 

TEE  Toxic  load  exponent 
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WWl  World  War  One 

Descriptors 
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